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Abstract 

Navigation of cells to the optimal environmental condition is critical for their survival and growth. Escherichia coli cells, for 
example, can detect various chemicals and move up or down those chemical gradients (i.e., chemotaxis). Using the same 
signaling machinery, they can also sense other external factors such as pH and temperature and navigate from both sides 
toward some intermediate levels of those stimuli. This mode of precision sensing is more sophisticated than the 
(unidirectional) chemotaxis strategy and requires distinctive molecular mechanisms to encode and track the preferred 
external conditions. To systematically study these different bacterial taxis behaviors, we develop a continuum model that 
incorporates microscopic signaling events in single cells into macroscopic population dynamics. A simple theoretical result 
is obtained for the steady state cell distribution in general. In particular, we find the cell distribution is controlled by the 
intracellular sensory dynamics as well as the dependence of the cells' speed on external factors. The model is verified by 
available experimental data in various taxis behaviors (including bacterial chemotaxis, pH taxis, and thermotaxis), and it also 
leads to predictions that can be tested by future experiments. Our analysis help reveal the key conditions/mechanisms for 
bacterial precision-sensing behaviors and directly connects the cellular taxis performances with the underlying molecular 
parameters. It provides a unified framework to study bacterial navigation in complex environments with chemical and non- 
chemical stimuli. 
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Introduction 

Living systems detect changes in the environment and try to find 
optimal conditions for their survival and growth. As one of the 
best-studied systems in biology, bacterial chemotaxis allows 
bacteria (such as Escherichia coli) to sense chemical gradients and 
navigate toward attractant or away from repellent [1—4]. This 
gradient sensing strategy makes cells move unidirectionally toward 
the extreme levels of stimuli. However, for other natural factors 
(such as pH and temperature), the physiological optimum does not 
locate at the extreme but at some intermediate level in the 
respective gradient. To find such intermediate point, it requires a 
more sophisticated strategy, namely, precision sensing. Both pH taxis 
[5-9] and thermotaxis of E. coli [10-16] provide us inspiring 
examples of precision sensing. 

Amazingly, E. coli cells use the same signaling system to achieve 
these different navigation tasks. Different external signals are 
sensed by several types of transmembrane chemoreceptors, among 
which the Tar and Tsr receptors are the most abundant [17]. For 
chemotaxis, binding of attractant (or repellent) molecules to 
chemoreceptors triggers their conformational changes and affects 
the autophosphorylation of the histidine kinase CheA [3,4]. 
Analogous to ligand binding in chemotaxis, both temperature and 
pH affect the conformational state of chemoreceptors and hence 
the CheA activity. Regardless of the way of being activated, 
phosphorylated CheA transfers its phosphate group to the 



response regulator CheY in the cytoplasm. The phosphorylated 
CheY molecules (denoted as CheY-P) then bind to the flagellar 
motors, increase their probability of clockwise rotations, and cause 
E. coli to tumble. The resulted alternating run and tumble pattern 
can steer cells to advantageous locations. To make temporal 
comparisons of stimuli, a short-term memory (or adaptation 
mechanism) is required [4,18,19]. This is achieved by the slow 
methylation-demethylation kinetics, as catalyzed by two enzymes 
(CheR and CheB) that add and remove methyl group at specific 
sites of receptors, respectively. 

How does a bacterium navigate through its environment with 
different chemical and nonchemical cues by using the same 
signaling and motility machinery? How do bacterial cells make 
decisions under competing chemical and/ or nonchemical signals? 
How accurately and reliably can bacteria find their favored 
conditions (such as the preferred temperature or pH)? and how do 
they tune their preference for precision sensing? We aim to address 
these questions under a unified theoretical framework, given that 
different taxis behaviors are based on the same sensory/motion 
machinery. 

To this end, we develop a multi-scale model which incorporates 
intracellular signaling pathways into bacterial population dynam- 
ics. The continuum population model reveals a simple theoretical 
result for the steady state cell distribution, which is found to be 
determined by the direction-dependent tumbling rates (transmit- 
ted through intracellular signaling pathways) as well as the 
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Author Summary 

Bacteria, such as E. coli, live in a complex environment with 
varying chemical and/or non-chemical stimuli. They con- 
stantly seek for and migrate to optimal environmental 
conditions. A well-known example is E. coli chemotaxis 
which direct cell movements up or down chemical 
gradients. Using the same machinery, E coli can also 
respond to non-chemical factors (e.g., pH and tempera- 
ture) and navigate toward certain intermediate, optimal 
levels of those stimuli. Such taxis behaviors are more 
sophisticated and require distinctive sensing mechanisms. 
In this paper, we develop a unified model for different 
bacterial taxis strategies. This multiscale model incorpo- 
rates intracellular signaling pathways into population 
dynamics and leads to a simple theoretical result regarding 
the steady-state population distribution. Our model can be 
applied to reveal the key mechanisms for different taxis 
behaviors and quantitatively account for various experi- 
mental data. New predictions can be made within this new 
model framework to direct future experiments. 



Here, we do not distinguish the methylation dynamics for different 
types of receptors (which appears to be regulated by the receptor- 
specific activity [30]) and consider m as the average methylation 
level of the whole receptor cluster. This treatment does not affect 
our main results since we are only interested in the total receptor- 
kinase activity, a, of the entire receptor cluster. 

A swimming bacterial cell may change its direction due to two 
mechanisms: the active transition to the tumbling state and the 
passive rotational diffusion (characterized by the rotational 
diffusion rate, Zq). According to the measured flagellar CW bias 
[31], the (instantaneous) rate of going into the tumbling state can 
be described as: z\(a) = x~ l {a/ Ki/ 2 ) H , where t^0.2s is the 
duration time of the tumbling state, K\ j2 represents the activity 
level at which the CW bias is 0.5, and //^10 denotes the 
ultrasensitivity of the motor response to CheY-P. Combining these 
two effects, the effective tumbling rate z(a) is given as: 

1 / \ H 

z(a) = z 0 + z l (a) = z 0 +-(^^J . (2) 



dependence of the swimming speed on external factors (such as 
temperature). This new finding enables us to systematically 
analyze bacterial navigation in chemical, pH, and temperature 
gradients. From each application, we have made quantitative 
comparison with the available experimental data and have gained 
new insights about the mechanisms of bacterial taxis. Our general 
model can be extended to study bacterial migration in complex 
environments (e.g., with a mixture of chemical and nonchemical 
stimuli) and provide quantitative predictions to be tested by future 
population level experiments. 

Results 

The Unified Model for Bacterial Taxis: A Simple Result for 
Steady State Cell Distribution 

Our unified model for bacterial taxis is developed on the basis of 
a number of previous models at different scales [20-29]. We 
incorporate microscopic pathway dynamics into the macroscopic 
transport equations [22,28] and derive a closed-form solution for 
the steady state cell distribution in chemical and/ or nonchemical 
gradients. In the following, we outline the main steps in obtaining 
this key result (Eq. 5), with more details about our model given in 
Text SI. 

The architecture of our model is illustrated in Fig. 1. The 
environmental signals (such as chemoattractants, pH, and 
temperature), denoted by S(x) = {S\,S2 ,.••}, can be sensed by 
different types of transmembrane chemoreceptors and converted 
into the total receptor-kinase activity, denoted by a. This activity 
represents the internal state of the cell and is described by the 
Monod-Wyman-Changeux (MWC) two-state model [23-26,29]: 



a(m,S) = 



1 



l + exp[7V/ a (m,S)]' 



(1) 



where N measures the degree of receptor cooperativity and f a 
represents the free energy difference between the active and 
inactive receptor conformations. The total activity, a, depends on 
the average methylation level of receptors, m, which restores a to 
the adapted level, do, over a time scale T a . For simplicity, the 
methylation rate is taken to be linear in a and hence the 

dm 

methylation dynamics can be described by: —— =(ao — a)/T a . 



In response to environmental signals, a population of bacteria 
will move in the physical space. Different from purely passive 
Brownian particles, cells also "distribute" in the internal state 
space, as each cell carries its own internal activity a when moving 
around. In the one-dimensional setup, let P-(x,a,t) denote the 
probability to find a cell being in the internal state a and moving in 
the " + " or " — " direction at (x,t). One can write down the 
master equation that governs the evolution of these probabilities 
(Text SI). As in many experiments, here we study the distribution 
of cells constrained in a finite chamber with a chemical or non- 
chemical gradient. The cell population distribution will equilibrate 
given enough time as the diffusion of cells balance the taxis drift 
effect. Using the zero-flux condition in the master equation leads 
to an exact expression for the steady state cell distribution p(pc) 
(Text SI): 



p(x) = 



Q 

(P + +P~)da=— exp 
v(x) 



(x')-z-(x') 
2v(x') 



dx' 



,(3) 
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Figure 1. Illustration of the multiscale architecture of our 
unified model. Environmental signals are sensed by the transmem- 
brane receptor-kinase complexes which controls the level of the 
intracellular response regulator (CheY-P). The response regulator 
controls the rotational direction of flagellar motors and the bacterial 
tumble frequency. Some environmental factors (such as temperature) 
can also affect bacterial swimming speed and motor switching 
dynamics. The population distribution of cells is finally shaped by the 
alternating tumble and swim behaviors. 
doi:1 0.1 371/journal.pcbi.l 003672.g001 
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where Q is the normalization constant and z-(x) = 



\zP ± da 



\P±da 

represents the average tumbling rate for the right or left moving 
cells at the same position x. 

It is clear from Eq. (3) that the cell distribution is determined by 
the two motility characteristics, tumbling rate (z) and swimming 
speed (v). On one hand, the local cell density is inversely 
proportional to the local swimming speed which may depend on 
the external condition S. Intuitively, it is easier for cells to leave a 
region if cells move faster there and thus cells spend more time in 
regions of low swimming speed. On the other hand, the cell 
density also depends on the accumulative (integrated) effect of the 
tumbling rate difference between the left and the right moving 
cells. For example, if z + (x) <z~(x), cells tend to move in the right 
( + ) direction on average because it is more difficult for cells to 
enter a region where they tumble more frequently. 

What is the origin of different tumbling rates (z + (x) ¥"Z~ (x)) for 
cells moving in different directions? The tumbling frequency is 
controlled by the CheY-P level which is proportional to the total 
activity, a. At any location x, the internal activity a is not fixed but 
distributed around its average, a. In fact, the average activity of the 
f aP + da 

left moving cells (a + = j ^ + — ) is different from that of the right 

_ \aP~da . 
moving cells [a = j ^ — ) as these two populations carry 

different average receptor methylation levels (memories). This 
activity difference can be evaluated (see Text S 1 and Figure S 1 for 

a "i \ a + - 2v dS (da\ 
more details : Aa = a —a « -^^—r~ -rzz , which is pro- 

portional to the (average) run length 2v/z and is valid as long as 
the adaptation time x a is much longer than the average run time 
l/z. The activity difference Aa can be used to evaluate the 
tumbling rate difference (Az = z + —z~)\ 



Azt 



Aa^lv 



dS /dlnz da 
dx V da dS 



(4) 



These functions can be determined by molecular and cellular level 
experiments, such as Ref. [26] for f a and Ref. [31] for z(a). Here, 
our model shows how population level (macroscopic) behaviors of 
cells can be predicted quantitatively based on these molecular level 
(microscopic) signaling and response characteristics. The general 
expression for steady state cell distribution, Eq. (5), provides a 
unified framework to systematically study diverse bacterial 
navigation behaviors in response to different chemical and non- 
chemical gradients, as will be shown in the following. We will also 
compare our theoretical results with the available experiments and 
make quantitative predictions on population-level behaviors of 
bacteria in more complex environments. 

Population-Level Sensitivity of Bacterial Chemotaxis 

We first apply our unified model to the case of bacterial 
chemotaxis. For simplicity, the chemical gradient (e.g. aspartate), 
denoted by is specifically sensed by one type of receptors 

(e.g. Tar) whose average activity can be described by the two-state 
model, Eq. (1). The free-energy difference between the active and 
inactive states is given by [25,29]: 



f a (m, [L]) =f m +/l = -E m (m- m 0 ) + In 



1 + 

l + [L]/K A 



(6) 



where f m (pi) and /lQL]) denote the methylation- and ligand- 
dependent contributions, respectively. The prefactor E m is the 
free-energy change per added methyl group, mo is a reference 
methylation level, and K A and Kj represent the dissociation 
constants for the active and inactive conformations, respectively. 

E. coli swimming speed (v) does not vary with the aspartate 
concentrations and is treated as constant here. By Eq. (5), it is easy 
to derive the steady-state cell distribution: 



p(x)«Qexp[ I7 («)/ L ([L](x))], 



(7) 



with a dimensionless factor rj that represents the effective 
sensitivity of the bacterial population to the environment: 



By using the above expression for (z + — z~) in Eq. (3), we 
finally obtain a simple expression for the steady state cell 
distribution p(x): 



f](a) = HN(\-a) 



z 0 + T-\a/K l/2 ) H 



(8) 



dlnz da 



dS 



(5) 



The equation for a is given in Text SI. 

This general expression (Eq. (5)) for the cell distribution is the 
main theoretical result of our paper. It shows that the steady state 
cell distribution is determined by two separable effects, the local 
effect of swimming speed v(S) and the accumulative effect of the 
gradient-dependent tumbling rates governed by internal signaling 
dynamics. In a previous work [20,21], a simple relation, 
p(x) = Q/v(x), was derived by assuming that the tumbling rate 
directly depends on the local environment factor. This treatment, 
however, did not take into account the cell's internal state or 
memory. Therefore, although the l/v dependence in Ref. [21] 
agrees with our Eq. (5), the integrated effect of the intracellular 
signaling dynamics was not identified or captured before. 

The intracellular signaling response to specific stimuli and the 
motor response to the response regulator are characterized by the 
free energy function f a (m,S), and tumbling rate z(a), respectively. 



The effective sensitivity rj is proportional to the signal 
amplification factors at both the receptor and the motor levels 
(i.e., H and N). It is dampened by the rotational diffusion (zq), 
because random collision of cells with the medium reduces 
directed chemotaxic motion. The dependence of rj(a) on the 
average activity a is nontrivial: on one hand, an increase of a could 
significantly boost the intrinsic tumbling [z\ oc a H ) and thus 
suppress the negative role of rotational diffusion; on the other 
hand, chemoreceptors become less responsive at a higher activity 
level a. Consequently, a has to be in a narrow optimal range in 
order to achieve high sensitivity. 

Our model allows for quantitative comparison with experi- 
ments. As shown in Fig. 2, the same functional form given in Eq. 
(7) can be used to fit the cell distribution data [32] in different 
attractant gradients. The coefficient rj inferred from experiments 



d\L\ 
dx ' 



(Fig. 2) appears to decrease with the gradient steepness Ql = 

indicating a higher population-level sensitivity in shallower 
chemical gradients. This trend agrees with the observation in 
our model that the average receptor activity a tends to increase 
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Figure 2. The cell density distributions for bacterial chemo- 
taxis. Different density profiles correspond to different gradients of 
varying steepness G L . Symbols represent the experimental data from 
Ref. [32], whereas lines denote the fitting of our Eq. (7) to the data. We 
used K A = 7>mM, ^ 7 = 18^M for MeAsp here. 
doi:1 0.1 371 /journal.pcbi.1 003672.g002 

with the gradient steepness in closed geometry. This increase in a 
is caused by the back flow of cells due to diffusion as p(x) is peaked 
at the boundary with the higher attractant concentration (see 
Fig. 2). Note that this is different from the open geometry case 
where the cell density is constant and there is no net diffusive back 
flow of cells. 

Tunability vs. Accuracy of E. coli pH Taxis 

E. coli can sense pH changes in the environment. According to 
recent experiments [8], Tar receptors exhibit an attractant 
response to a decrease of pH while an opposite response was 
observed for Tsr. The balance between the two opposing receptors 
leads to a preferred pH level for the wild-type E. coli, i.e., precision 
sensing, as suggested by our recent model study [9] of intracellular 
pH responses. Here, we examine how accurately and how robustly 
a population of bacteria could find their preferred pH level. 

The extracellular pH modulates the receptor-kinase activity 
primarily by affecting the periplasmic domain of the Tar and Tsr 
receptors. The total receptor-kinase activity can be described by 
the generalized MWC model for heterogeneous types of receptors. 
The total free energy between the active and inactive states is 
given by 



fa(m,pH) =f m + rxfi (pH) + r^CpH) 

l + 10^" pH 



~fm 



0=1,2 1- 



(9) 



here. Using Eq. (5), one can easily obtain the cell distribution in a 
pH gradient: p(x)^Qexp[— Feff(x)], with the effective potential 
Feff = — ^ [ri/i (pH) + 7*2/2 (pH)]. The competition between the 
opposing pH dependence of f\ (from Tar) and fi (from Tsr) leads 
to the accumulation of cells at an intermediate (preferred) pH 
level, which can be analytically determined by the condition 
Feff^pH*) = 0; see Text SI for more details. 

The preferred pH is mostly sensitive to the relative abundances 
of receptors {ri/fy) and the values of Kf and Kf. Based on our 
analysis and simulations, we find an empirical equation (Text SI), 



pH* ; 



— -^logio, ,. 



(10) 



which shows the logarithmic dependence of the preferred pH on 
the relative abundance of Tar and Tsr (Fig. 3A). The coefficient X 
in Eq. (10) varies with the dissociation constants and can be 
interpreted as a measure of tunability of the preferred pH upon 
changing the Tar/Tsr ratio. Theoretically, this coefficient is close 
to one when Kf&Kf (Text SI). Numerically, we also 

found that higher tunability (A>1) can be achieved if Kf <Kf 
and the opposite holds for Kf >Kf (Fig. 3A). Our results can be 
compared with the recent experiment [8], where the pH 
preference point (pH*) was observed to shift from 8.0 to 7.5 when 
the Tar/Tsr ratio changed from 0.43 to 1.25 (symbols in 

Fig. 3A). Using these data with the empirical Eq. (10) yields 
(Kf + Kf)/2 = 1.6 and X =1.08, which together indicate that 
Kf^Kf^l.6. 

In addition to the preferred pH, our population model also tell 
us the dispersion (or accuracy) of bacteria seeking and aggregating 
around their favored pH, which can be quantified by the standard 
deviation of the cell distribution in the pH scale. As shown in 
Fig. 3B, the dispersion measure turns out to be minimal for the 
scheme Kf = Kf , compared to the dispersion for either Kf > Kf 
or Kf <Kf . Therefore, one possible advantage for having 
Kf ~Kf (Fig. 3A) in E. coli is the optimal accuracy of pH sensing 
at population level, though with a tradeoff of a modest tunability 



: 1) for Kf 1 



■10"? 




B, 



-8 0.21 



2 0.19 



0.17 



0 0.5 
.(Tar/Tsr) 




-0.2 0 0.2 
log (Tar/Tsr) 



where r\ and ri denote the fractions of Tar (q= 1) and Tsr (q = 2) 
in the receptor cluster, respectively. The dissociation constants 
K[ 2 an d Kf 2 for the inactive and active receptors are expressed in 
the pH scale. The observed opposite responses to pH changes 
indicate that Kf <K[ for Tar and Kf >K[ for Tsr. Without loss 



= 9.0 and 



6.0 for numerical 



of generality, we set 
examples. 

As the E. coli motor speed does not vary significantly with the 
external pH [33], we can take the swimming speed as constant 



Figure 3. Tunability and accuracy of bacterial pH taxis. (A) The 

preferred pH* versus the logarithm of the Tar/Tsr ratio to base 10 for 
three representative parameter regimes: Kf=Kf, Kf<Kf, and 
Kf>Kf. The red symbols represent the experimental data from Ref. 
[8] and seem to coincide with the model curve for Kf=Kf. (B) The 
standard deviations of the cell distributions as a function of log 10 (r\/r 2 ) 
for the three representative parameter regimes: Kf =Kf, Kf <Kf, and 
Kf>Kf. In the above numerical examples, we have fixed 
(Kf+Kf) 12 = 1.6, K[ = 9.0 and A* = 6.0. 
doi:1 0.1 371/journal.pcbi.l 003672.g003 
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The Inversion of Thermal Response at the Critical 
Temperature 

Bacteria are able to sense thermal variations and migrate 
toward their favored temperature [10-15], another example of 
precision sensing. However, unlike in pH sensing where two types 
of receptors, Tar and Tsr, respond in opposite ways to a pH 
change, temperature sensing can be achieved by a given type of 
receptor (Tar) which changes the sensing mode (from being a 
warm sensor to a cold sensor) as its methylation level m increases 
across a critical level m c (&2) [12,13]. Added to the complexity is 
the fact that temperature affects many other aspects of motility, 
such as the swimming speed and the motor switch sensitivity. 
Here, we first demonstrate how a chemoreceptor acts as a thermal 
sensor that inverts response at some critical temperature. In the 
next section, we will study how all those temperature-sensitive 
factors affect thermotaxis at the population level. 

For simplicity, we consider E. coli cells that only express Tar 
receptors and migrate in a linear temperature gradient. In general, 
the total free energy for the Tar activity can be described as 



f a (m,T)~ -E m (m-m 0 )-(m-m c yg(T), 



(11) 



where g(T) describes how temperature affects the total free energy 
and where E m refers to the free energy change per added methyl 
group (in units of kT) at a given reference temperature Tq (i.e., 
g(To) = 0). Note that a linear function g{T) = a{T — Tq) with a > 0 
was used in a previous model of thermotaxis [16]. However, it is 
easy to verify that as long as g'(T)>0, the Tar receptor switch 



A 



£ -0.02 
P 

"g -0.04 



warm sensor 


cold sensor 


/ T c 

, , — 


— m(T) - m c 

- - da/dT 


■ 



36 38 
T(°C) 




da 

from being a warm sensor (-^ < 0) when m<m c to a cold sensor 
da 

(— >0) when m>m c ; see Fig. 4A. 

It is observed experimentally that the adapted activity also 
changes with temperature [26,34]. This can be modeled by 
a 0 = 1/ [1 . For E. coli, it is reported that a 0 « 1/3 at 

room temperature T = 22°C and ao « 1/2 at To = 32° C, leading to 
an estimate of /?~0.07/°C. The critical temperature T c at which 
Tar inverts its response is defined by m(T c ) = m c . Using the 
condition a = ao, we can obtain the steady state methylation level 
m(T) as well as the critical temperature: 
T c = Tq + NE m (m c — mo) / which is determined by the (up- 
stream) receptor kinetics. This also leads to a simple relationship: 



m(T) — m c 



P(T-T C ) 
N[g(T)+E m 



(12) 



showing that the average methylation level relative to m c changes 
sign at the critical temperature T c (Fig. 4A). According to Eq. (12), 
when T < T c , the Tar methylation level is less than m c (m<m c ), so 
the Tar receptor is a warm sensor driving the cell towards T c from 
lower T; when T>T C , the Tar methylation level is greater than 
m c (m>m c ), now the Tar is a cold sensor driving the cell towards 
T c from higher T. This is the basic mechanism for cell 
accumulation around T=T C . 

Two Channels Drive Bacterial Thermotaxis: Speed and 
Sensing 

Besides the receptor-kinase activities, temperature also affects 
other aspects of the system. For example, it is observed for E. coli 
that both the motor dissociation constant, K\/2, and the swimming 
speed, v, change with temperature [34,35]. It remains unclear 
whether and how different temperature-sensitive factors affect the 
performance of bacterial thermotaxis. 

Using Eq. (5), one can derive the steady-state cell distribution 
over the temperature range [T_,T + ]: 



Q 



Q 



v(T)Z(T) v(D 



exp 



N 



(T-T c )rj(f)6(f)df 



,(13) 



Figure 4. Inverted response to temperature changes and 
bacterial thermotaxis. (A) The steady-state methylation level 
subtract the critical methylation level, m(T) — m c , and the receptor 
response to temperature changes, da/dT, as a function of temperature. 
The critical temperature T c is determined by the crossing point where 
m(T c ) = m c (or equivalently da/dT = 0). Tar acts as a warm sensor for 
T <T C and a cold sensor for T>T CI which drives the cells towards T c 
from both sides. (B) The steady-state cell distribution, p, as a function of 
temperature. For illustrative purposes, we assume that the swimming 
speed, v, increases linearly with temperature and that the motor 
dissociation constant is K l/2 (T) = l/[l +e~ y(T ~ To) ], with a constant 
parameter y. Three cases are considered. The red solid line corresponds 
the case where both v and K\/ 2 are constant; the blue dot-dashed lines 
is for the case of constant K\/ 2 (i.e. y = 0) and dv/dT=\.0; the green 
dashed line is generated by using y = 0.1/°C and dv/dT = 1.0. Evidently, 
the steady-state cell distribution can be changed by the temperature 
dependence of the speed v(T), but it is insensitive to the temperature 
dependence of motor sensitivity K\/ 2 (T). Here, we fix 9 = 0.1 in all 
numerical examples. Other parameters used include: 9 = 0. 1/°C, 
£ = 0.07/°C, E m = \.0, H=\0, N = 6, m 0 = 1.94, m c = 2.0, T 0 = 32°C, 
z 0 = 0.28rad 2 /s, and t = 0.2s. 
doi:1 0.1 371 /journal.pcbi.1 003672.g004 



where the function Z{T) represents the effect of the direction- 
dependent tumbling rates governed by the thermosensory system. 
In Eq. (13), f](T) is the effective sensitivity defined in Eq. (8) and 
has weak dependence on temperature through both ao(T) and 

g\T) 

K]n(T). The function 0(T) = —— is introduced for 

1 ' g(T) + E m 

convenience and represents the effect of the sensory system. 

Eq. (13) shows that there are two independent channels affecting 
bacterial thermotaxis: one is the swimming speed v(T), and the 
other is the sensory system that controls the rotational direction of 
flagellar motors. In contrast to the local speed effect v(T) which is 
direct and memoryless, the sensing effect Z(T) is indirect 
(channeled through signaling networks and motor control) and 
relies on the slow adaptation dynamics which encodes memory for 
the system to sense the environment [19,29]. 

Near the critical temperature T c , we can compute Z(T) by 
keeping only the leading order term in (T—T c ), that is, 
Z(T)k Qxp[^(T c )6(T c )(T-T c ) 2 /(2N)}. The expression for the 
cell density follows: 
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P{T)r 



n 

AT) 



exp[-C(T-T c ) 2 ], 



(14) 



where C = 



Prj(T c )6(T c ) 
IN 



is a positive constant when /?>0 and £ 



6(T C )>0. It is clear from the above equation that for a constant 
v( T) cells will accumulate around the temperature T* = T c 
(Fig. 4B). The accumulation temperature T* can be shifted from 
T c by the dependence of the swimming speed on temperature, 
e.g., if v increases with temperature, cells tend to spend more time 
in regions of lower speed and thus aggregate at a lower 
temperature, i.e. T* <T C (Fig. 4B). 

The shift of T* from T c only weakly depends on K\/2, which 
indirectly affect C in Eq. (14) through r] (see Eq. (8)). In fact, even 
the shape of the distribution p(T) is not sensitive to Ki^(T), as 
shown in Fig. 4B. The insensitivity of thermotaxis to K\/2(T) is 
due to the fact that K\/2 only depends on the local temperature 
and is the same for different cells at a given position x, regardless 
of their direction of motion. In other words, K\/2(T) does not 
contribute to the tumbling rate difference Az(x) that drives the 
directed migration (taxis) of cells. 

Thermotaxis in Shallow Temperature Gradients: Model 
vs. Data 

Our theory can help explain some recent experiments 
measuring cellular behaviors in shallow temperature gradients 
[36,37]. It was observed that even the mutant bacteria lacking 
all chemoreceptors are still able to migrate toward high 
temperature [36], showing that there is an additional channel 
(other than sensing) in regulating bacterial thermotaxis. When 
the sensing channel (i.e. the bacterial signaling machinery 
which translates temperature stimuli into tumbling bias) does 
not work (e.g., due to deletion of the receptors), the 
temperature-dependent swimming speed can still cause the 
directed cell migration [36]. This is consistent with earlier 
work [20,21] and our model where mutant strains lacking all 
receptors can be described by 6(T) = 0 which leads to 
p&Q/v(T) as described by Eq. (13). 

In Ref. [37], the swimming speed v(T) for wild-type E. coli (with 
functional chemoreceptors) was measured at different tempera- 
tures. The speed profile appears to be a quadratic function of 
temperature and reaches its maximum at T = 30°C. We have 
quantitatively compared the inverse speed profile, p~Q/v(T), 
with the cell density data in Ref. [37] and found that the inverse 
speed profile alone could not account for the observation, 
especially the significant aggregation of cells at high temperatures 
(Fig. 5). This suggests that the thermosensory system may not be 
completely silent in shallow temperature gradients as suggested in 
[37]. We test this hypothesis by using Eq. (13) with the measured 
v(T) and the assumption 0^O.2/°C for T>35C. It turns out that 
our model, which includes both channels (speed and sensing), 
provides a better agreement with the observed data (Fig. 5). This 
result suggests that the thermosensory system may be active even 
in shallow temperature gradients. Further experiments are needed 
to examine and quantify the interplay between these two channels 
(speed and sensing) in shallow temperature gradients. 

Navigation under Opposing Chemical and Thermal 
Gradients 

Our unified model can be applied to study and predict bacterial 
taxis behaviors in more complex environments. As a final 
example, we investigate the behavioral response of the Tar-only 
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Figure 5. Comparison between model results and experimental 
data for £ coli thermotaxis. The blue symbols represent the cell 
density data obtained in Ref. [37] at ~42 min after applying a shallow 
temperature gradient (~0.02°C//mi). The black dashed line is the 
inverse speed profile, p(T) = Q/v(T) where v(T)= -0.09(r-30) 2 + 24 
is a quadratic fitting to the measured swimming speed in Ref. [37]. The 
red solid line corresponds to model Eq. (13) with 0 = O.1[1 + 
tanh(r-35)] and v(T)= -0.09(r-30) 2 + 24. Other parameters used 
are the same to those in Fig. 4B. 
doi:1 0.1 371/journal.pcbi.l 003672.g005 



mutant cells over the interval [xo,Xi] under a temperature gradient 
( —j- = Qt > 0) and an opposing chemoattractant gradient 

(— U = — (?l<0). The total free energy is modified by adding 

an additional ligand-depend free energy /l(x) to Eq. (11). In this 
case, the steady-state cell density is found to be (Text SI): 



p(x)? 



Q 

vZ 



exp 



df L 



(15) 



where the term Ql T ^ n describes the chemotactic drift, and the 
d[L\ 

other term GtO/l captures the interaction between the chemical 
and thermal signals. 

In the absence of attractant (i.e., i /l = 0), Eq. (15) recovers 
Eq. (13) for thermotaxis. Interestingly, for a uniform chemical 
background (i.e.,/^>0 and dfL/dx = 0), the interference effect 
QT0(T)fL([L])>0 tends to suppress the accumulation of cells at 
high temperatures. Quantitatively, a uniform chemoattractant 
background can shift the preferred temperature from T* = T c 
to a lower temperature T* = T c — N/l/ fi- When there is an 
attractant gradient (so that df^/ 'dx<0) imposed against the 
temperature gradient, the accumulation point can be shifted 
further. Specifically, as the chemical gradient steepens (Ql 
increases), the chemotacic response of bacteria become 
stronger, leading to a positional shift of their aggregation 
toward lower temperatures, as shown in Fig. 6. This example 
demonstrates the general capability of our model in making 
quantitative predictions on bacterial behaviors in complex 
environments (with multiple and competing chemical and 
nonchemical stimuli), which can be used to guide future 
experiments. 
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0.15 




Figure 6. The cell density profiles under two opposing 
chemical and thermal gradients. The temperature gradient used 
here is from 12 C to 42° C in a channel of length 1.0cm. Different 
density profiles correspond to different attractant (MesAsp) gradients 
Gx = 0.0, 3.0, 5.0, and 6.0 /jM/cm but the same concentration at the 
middle point: [L](x = 0.5cm) = 3.0^M. Other parameters used are the 
same to those in Fig. 4B. 
doi:1 0.1 371 /journal.pcbi.1 003672.g006 

Discussion 

In this paper, we have incorporated the intracellular signaling 
pathways into the bacterial population dynamics and developed a 
unified model to study bacterial navigation in chemical and 
nonchemical gradients. This model leads to a general result, which 
shows that the steady state cell density is determined by the 
accumulative effect of the direction-dependent tumbling rates as 
well as the local swimming speed. In the following, we discuss 
some of the specific findings and related possible future directions. 

The Push-Pull Mechanism for Precision Sensing 

From the population model, we can construct an effective 
potential function V e ff oc — In (p), which provides a useful scheme 
to visualize different cases of bacterial taxis, as summarized in 
Fig. 7. The effective potential for chemotaxis decreases monoton- 
ically with the chemoattractant concentration and thus steer cells 
up the chemical gradient (Fig. 7). Our application to pH taxis 
illustrates how the competition between two pH sensors (Tar and 
Tsr) determines the preferred pH for the wild-type cells expressing 
both Tar and Tsr: a push-pull mechanism here creates a potential 
well for bacteria to accumulate (Fig. 7). In the case of E. coli 
thermotaxis, the push-pull mechanism is more subtle as the 
"push" and the "pull" are provided by the two sub-populations of 
Tar receptors with their methylation levels above or below the 
critical level (m c ). This leads to a well-defined critical temperature 
where cells tend to accumulate (Fig. 7). The push-pull mechanism 
is likely a general strategy for precision sensing. For example, it 
was found that two receptors, Tar and Aer, leads to a preferred 
level of oxygen for E. coli aerotaxis [38], which may also be studied 
within our unified model. 

Robustness and Sensitivity 

E. coli chemotaxis has served as a model system in studying 
robustness of biochemical networks [34,39,40]. Bacteria exhibit 
thermal robustness in their chemotaxis network output by 



chemotaxis thermotaxis 




pH T 



Figure 7. Schematic illustration of the effective potential 
for chemotaxis, pH taxis, and thermotaxis. In the case of 
chemotaxis, F e ff([L]) decreases monotonically as the chemoattractant 
concentration [L] increases. For pH taxis, K e ff(pH) decreases with pH 
for Tsr-only mutant cells and increase with pH for Tar-only mutant. 
Based on the push-pull mechanism, V e ff(T) for the wild-type E. coli 
represents the balancing effect between Tar and Tsr, leading to a local 
minimum in the effective potential. In the case of thermotaxis, V&{T) 
can be shifted by the effect of temperature-dependent swimming 
speed v(T). It is, however, insensitive to other temperature effects 
such as the temperature dependence of motor response, Ki/ 2 (T) = 
l/[l -\-e~ y ( T ~ To) ], parameterized by y. 
doi:1 0.1 371/journal.pcbi.l 003672.g007 



counterbalancing temperature effects on different opposing 
network components [34]. For example, the dissociation constant 
K\/2(T) for the motor switch is observed to increase with 
temperature [35] . This effect balances the increase of the adapted 
CheY-P level with temperature such that the motor switch is able 
to operate in a narrow optimal range with ultrasensitivity [31]. 
This, however, raises a question for bacterial thermotaxis: do those 
temperature-sensitive factors, such as K\/2(T) and ao(T), hinder 
the thermotactic performance? According to our model analysis, 
the steady state distribution of cells in a temperature gradient is 
mainly determined by two effects: the temperature-dependent 
swimming speed and the direction-dependent tumbling rates. The 
system is actually robust/insensitive to those instantaneous/local 
temperature-sensitive factors (e.g. K\/2(T)) which do not contrib- 
ute to the tumbling rate difference at any spatial point. The 
insensitivity of thermotaxis to Ki^T), as shown from our model, 
is a highly desirable feature of the system as it allows robust 
thermotaxis without sacrificing motor-level sensitivity. 

Navigation in Complex Environments 

In the natural environment, cells are often exposed to multiple 
chemical stimuli [41]. Our general model can be applied to study 
such cases (with a specific example discussed in Text SI). The 
density of cells subject to a multitude of chemical gradients shall be 
given by p(x) « Q exp (rj Y^ q r qfq), where f q ([L] l9 [L] 2 ,...) denotes 
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the free energy contribution from all the chemical signals that are 
sensed by the type-g receptors. Quantitative predictions can be 
made about how bacterial cells integrate and respond to mixed or 
competitive chemical signals and how their response changes with 
the composition and relative abundance of their sensors. More 
complex situations exist when different stimuli are interdependent 
and/or interfere with non-chemical factors. For example, the 
chemical environment can be modified through consumption and 
secretion by the bacteria, a dynamical process depending on the 
bacterial cell density [36] . In addition, temperature can change the 
metabolic rates of bacterial cells and create temperature-depen- 
dent chemical (nutrients, oxygen) gradients. How cells navigate 
under such complex circumstances and how such behaviors lead 
to survival/growth benefits remain unclear. Our model can be 
extended to study those phenomena and help address those 
fundamental questions. 

In sum, the work presented here provides a general model 
framework to study population behaviors in the presence of both 
chemical and non-chemical signals based on realistic intracellular 
signaling dynamics. 
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